Non-mean-field theory of anomalously large double-layer capacitance 
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Mean-field theories claim that the capacitance of the double-layer formed at a metal/ionic con- 
ductor interface cannot be larger than that of the Helmholtz capacitor, whose width is equal to the 
radius of an ion. However, in some experiments the apparent width of the double-layer capacitor 
is substantially smaller. We propose an alternate, non-mean-field theory of the ionic double-layer 
to explain such large capacitance values. Our theory allows for the binding of discrete ions to their 
image charges in the metal, which results in the formation of interface dipoles. We focus primarily 
on the case where only small cations are mobile and other ions form an oppositely-charged back- 
ground. In this case, at small temperature and zero applied voltage dipoles form a correlated liquid 
on both contacts. We show that at small voltages the capacitance of the double-layer is determined 
by the transfer of dipoles from one electrode to the other and is therefore limited only by the weak 
dipole-dipole repulsion between bound ions, so that the capacitance is very large. At large voltages 
the depletion of bound ions from one of the capacitor electrodes triggers a collapse of the capaci- 
tance to the much smaller mean-field value, as seen in experimental data. We test our analytical 
predictions with a Monte Carlo simulation and find good agreement. We further argue that our 
"one-component plasma" model should work well for strongly asymmetric ion liquids. We believe 
that this work also suggests an improved theory of pseudo-capacitance. 



I. INTRODUCTION 

The rising demand for compact forms of energy storage 
with high power output has resulted in increased interest 
in electrochemical capacitors (ECs) [l|, 0]. An EC is a 
pair of metal electrodes separated by an ionic conductor, 
such as an aqueous solution of ions, an ionic liquid JSj], a 
super- ionic crystal jj] , or an ion-conducting glass [5|, |6| . 
ECs with extremely high area per unit volume ("super- 
capacitors") already have a large number of applications. 
In this paper our focus is not on large surface area, but 
on the deeper physical question of a maximum possible 
capacitance per unit area. 

In a conventional double-plate capacitor, where metal 
electrodes of area S are separated by an insulator of 
width d and dielectric constant e, the capacitance C = 
eS/A'Kd (in Gaussian units). In an EC, the intervening 
medium is actually a conductor with finite conductivity 
CT, but with blocking of both ionic and electronic cur- 
rent at the electrode interface. The relation C = eS/iird 
is therefore only valid at sufficiently high frequencies 
w ^ Aira/e, where the bulk of the ionic medium behaves 
as an insulator. We concern ourself with the opposite 
limit io <^ Aira/e, where polarization of the ionic medium 
eliminates electric field in the bulk and the capacitance 
of the EC is determined by the formation of thin electro- 
static double-layers (EDLs) at both electrodes. 

How large can the capacitance be for these double- 
layers? The commonly-accepted expression for the max- 
imum possible capacitance of an EDL goes back to 
Helmholtz [7| , who assumed that the charge of the metal 
surface is compensated by a layer of counterions with di- 
ameter a residing on the surface of metal. The resulting 
"Helmholtz capacitance" is given by 



Ch = eS/2Tra. 



(1) 



For a double-plate capacitor, where the EDLs formed at 
both electrodes can be thought of as two equal capaci- 
tances connected in series, the maximum capacitance is 
Ch/2 = eS/i-Ka. For a = 2A and e = 5, as we use below, 
Ch/2S = 22 ^iF/cm^. 

A recent experiment [6|, however, has reported much 
larger values of the EDL capacitance in phosphosilicate 
glasses placed between platinum electrodes (see Fig. [S]). 
Capacitance per unit area as large as 400 /iF/cm was 
measured, corresponding to an effective capacitor thick- 
ness d* = eS/iirC in the range 0.2 ~ 0.7 A, much smaller 
than any ion radius. The glass was held at a tempera- 
ture of 573 K, at which only the smallest ions, Na+ with 
diameter a — 2A, are mobile. The dielectric constant of 
the glass e is between 2 and 10. 

Current theories of EDL capacitors, based on the 
mean-field approach, fail to explain such large capaci- 
tance. The most widely-used theory is that of Gouy, 
Chapman, and Stern (CCS) |8l-[l0|. which extends the 
Helmholz capacitor concept to allow for the thermal mo- 
tion of counterions. In this approach, neutralizing ionic 
charge is imagined as a stack of thin uniform layers placed 
parallel to the charged electrode, with the charge density 
of each layer dictated by the Poisson-Boltzmann equa- 
tion. When excluded volume of ions is properly taken 
into account |ll| - [la | such theories lead only to a larger 
effective capacitor thickness, and therefore a smaller ca- 
pacitance than the Helmholtz value. 

In this paper we propose an alternate theory to ex- 
plain the large differential capacitance of the EDL. We 
abandon the mean-field approximation and deal instead 
with discrete ion charges, which interact strongly with 
the metal surface in a way that is not captured by the 
mean-field approximation. For simplicity, this paper fo- 
cuses primarily on the case where only the cations are mo- 
bile, and therefore may form the EDL, while anions com- 



prise a fixed background of negative neutralizing charge. 
This "one component plasma" (OCP) model is used here 
to describe the ion-conducting glasses examined in Ref. 
[g . We also suggest a number of other systems to which 
it can be applied. 

One other application of the OCP model is to super- 
ionic crystals, where only the smallest positive ionic 
species (such as Na+ or Li+) is mobile. Less obvious is 
the application of the OCP model to ionic liquids, where 
both positive and negative ions are mobile. Nevertheless, 
we show below that for ionic liquids composed of monova- 
lent cations and much larger, rigid anions, one can think 
about space-filling anions as a weakly-compressible nega- 
tive background on which the small cations rearrange to 
form the EDL. If the cations are multivalent, then mono- 
valent anions can again be treated as a negative back- 
ground even if they are of the same size as the cations. 
These applications will be explained in greater detail in 
a later section. 

We begin our theory of the OCP model by noting that 
a cation adjacent to a metal electrode produces electronic 
polarization of the metal surface, and the cation experi- 
ences an attraction to the resulting image charge. For 
ions of small radius, the image attraction is significantly 
larger than the thermal energy fc^T, so that ions form 
stable, compact ion-image dipoles at the metal surface. 
Cations may also experience some chemically-specific at- 
traction to the metal electrodes, which enhances the ef- 
fects of image attraction. Repulsion between adjacent 
dipoles results in the formation of a strongly-correlated 
liquid of dipoles along the surface of both electrodes. 

The adsorbed cations leave behind a region of nega- 
tive background, or a depletion layer. Thus, each border 
of the sample is spontaneously polarized in the direc- 
tion perpendicular to it (see Fig. [T]). When a voltage V 
is applied between the electrodes, these opposite-facing, 
spontaneous polarizations are easily rearranged in the di- 
rection of the external electric field, leading to the rapid 
build-up of an electronic charge ±Q on the metal sur- 
faces. 

As the charge ±Q is added to the electrodes, cations 
are driven to detach from the positive electrode and to 
bind to the negative one. Below, we demonstrate that 
the resulting capacitance dQ/dV of the two EDLs can 
be significantly larger than the Helmholtz value Ch/2, 
since the dipole-dipole repulsion that resists ion transfer 
is relatively small. In other words, the effective capacitor 
thickness d* can be much smaller than the ion radius. Be- 
low we derive an expression for d* that is reasonably close 
to experimental values. Our theory also explains another 
peculiarity of the experiment [6| , namely the sharp drop 
in capacitance at a certain critical voltage (as in Fig. 
|6]). We show that, indeed, the capacitance should col- 
lapse to a much smaller value at a particular nonzero 
voltage. Contrary to the mean-field theories developed 
in Refs. |lll - [l6| |. this drop-off (shoulder) in capacitance 
is not driven by excluded volume effects among bound 
ions. Rather, it occurs far below the complete filling of 




FIG. 1: A capacitor consisting of parallel metal plates (lightly 
shaded) bounding an ionic conductor with mobile positive 
ions. The neutral region of the ionic conductor is heavily 
shaded, while negative depletion regions {zi and Z2) are left 
white. The relative size of the ions has been greatly exagger- 
ated. 



an ionic layer at either electrode [l7|, when the voltage 
difference induces the positively-charged electrode to lose 
all of its bound ions. 

The following section develops the OCP model of ca- 
pacitance for a parallel double-plate capacitor at low tem- 
peratures, where all entropic effects can be ignored. Sec- 
tion III extends this low-temperature theory to the case 
of an asymmetric double-plate capacitor, where the an- 
ode and cathode have different areas. In the limit of 
one area going to infinity, this includes the description 
of a single EDL capacitance. Section IV examines the 
temperature-dependence of the capacitance, and section 
V presents the results of a simple Monte Carlo simulation 
and compares them with our analytic theory. Section VI 
examines the application of the OCP model to ionic liq- 
uids, and section VII explores its relevance for electrolyte 
solutions in water. In Section VIII we summarize our 
main conclusions and we argue that our paper suggests 
an improved theory of "pseudo-capacitance" . 

A short version of this paper was published in Ref. [l^l . 



II. LOW-TEMPERATURE THEORY OF A 
SYMMETRIC PARALLEL-PLATE CAPACITOR 



We consider the case of a parallel-plate capacitor, 
where an ion-conducting medium is bounded on two sides 
by metal plates which are maintained at a relative poten- 
tial difference V by some voltage source. The interven- 
ing medium is modeled as a fixed negative background 
with charge density —eN, upon which resides a neutral- 
izing concentration of mobile positive ions with charge 
+e and bulk density N. We model the conducting ions 
as hard spheres with diameter a, and this paper gener- 
ally assumes that Na^ <C 1. When an ion is up against 
the metal surface, it experiences an attractive potential 



energy of approximately 



e 
2ea' 



(2) 



For T = 573 K, e = 5, and a = 2 A, as estimated for the 
experiments of Ref. §, we get \uirn\lkBT ~ 15, so that 
such ions are bound strongly to the surface. In fact, the 
attraction of ions to the metal surface may be stronger 
than Uim, either because of local behavior of the dielectric 
constant or because of some specific chemical affinity that 
ions may have for the metal surface. In general, we will 
write the strength of ion attraction to the metal surface 
as ^Uirm where 7 is a positive constant of order one. 
Unless otherwise stated, all numerical estimates will use 

7 = 1- , 

At a given voltage V , some area densities ni and 712 
of ions bind to the anode and cathode, respectively. The 
attachment of these positive ions to the metal surface re- 
sults in the formation of regions with negative net charge 
qi and 92 near the anode and cathode, respectively. Each 
of these charges exactly cancels the net charge of the ad- 
joining plate and its bound ions, so that there is no elec- 
tric field in the bulk of the ionic conductor. This implies 
qi = —eniS — Q and 52 = —en2S + Q, where Q is the 
electronic charge that moves through the voltage source. 
In other words, the dipoles at each metal-glass interface 
are effectively embraced by a capacitor composed of the 
charge q and its positive image — g in the metal. As in ev- 
ery plane capacitor, the charges gi and (72 are uniformly 
distributed along the plane. 

Since positive ions gain a large energy \uim\ by adsorb- 
ing to the metal surface, at equilibrium there must be a 
correspondingly large potential difference Vim = \uim\/e 
between the metal surface and the bulk of the ionic con- 
ductor, so that the chemical potential of ions is uniform. 
This large potential difference can be defined as the po- 
tential of zero charge (PZC), which has a nonzero value 
because of the strong spontaneous polarization of the 
system. In this paper we consistently use the external 
electrode-electrode potential V, which by definition is 
zero when there is no charge Q on the electrode surfaces, 
to define capacitance, rather than the potential difference 
between the electrode and the bulk. 

With such a large internal potential difference Vim ^ 
ksT/e at each electrode, the negative regions are 
strongly depleted of ions even at very small applied volt- 
age V. The charges qi and 52 therefore constitute de- 
pletion layers of width zi and Z2 which form at the an- 
ode and cathode, respectively; here it is assumed that 
2i, 02 ^ o- These layers are devoid of mobile ions and 
have a charge density equal to that of the negative back- 
ground, so that qi — —eNzi and (72 = —eNz2- Thus 



eNzi — eni = Q/S 
eNz2 — en2 = —Q/S, 



(3) 
(4) 



where Q is the electronic charge that moves through the 
voltage source of the capacitor. The electrostatic energy 



associated with the formation of the depletion layers and 
their corresponding positive image charge in the metal 
can be estimated as 



Ude 



2ne'^S 
3eN 



ni 



Se 



n2 



Se 



(5) 



Fig. [2] gives a schematic depiction of the potential en- 
ergy w{z) of an ion near the anode as a function of 
its distance z from the metal surface, relative to a po- 
sition in the bulk of the ionic conductor. There are 
two main contributions to w{z): the attraction energy 
Wim(z) = —e^/'iez of the ion to its image charge and the 
electrostatic energy Wdep{z) = 2ne'^N{zi — z)'^/e associ- 
ated with moving the ion from the bulk into the depletion 
region. At equilibrium, these two contributions satisfy 
Wmi(a/2) = Ui„i = — W(iep(a/2), so that bound ions at 
the surface have the same energy as ions in the bulk and 
the chemical potential is uniform. 




FIG. 2: Schematic depiction of the potential energy 'w{z) of 
an ion in the vicinity of the metal plate. The image attraction 
energy Wim and the energy associated with moving through 
the depletion layer Wdep are shown as a function of distance 
z from the metal surface, along with the total w = Wim + 
TOdep. At equilibrium, ions bound to the surface have the 
same potential energy as ions in the bulk. 



Eq. dni) assumes that positive bound ions are effec- 
tively neutralized by their negative image charges, so that 
these ions can be ignored in calculating the electrostatic 
self-energy of the depletion region and its image charge. 
One can consider, however, the small correction to Eq. 
([5]) resulting from the interaction between bound ion- 
image pairs and the negative background. In the limit 
zi,Z2 ^ a, this correction can be viewed as the poten- 
tial energy of a collection of ion-image dipoles, each with 
dipole moment \p\ = ea, aligned with an external field E 
produced by the negative depletion layer and its positive 
image charge. If the number of dipoles on a given metal 
surface is nS, then the total interaction energy of these 
dipoles with the external field is —nSp- E, so that to the 
energy associated with forming the depletion layers there 



should be added a term 



K 



dE 



-A-n-u.mS 



ni Til 



eS 



+ ^2 «2 



Q_ 

eS 



(6) 

In addition to the ions' interaction with the back- 
ground, there is a positive energy associated with repul- 
sion between bound ions. When the density of bound 
ions is low enough that na^ <^ 1, ions repel each other 
by a dipole-dipole interaction: the potential created by a 
bound ion and its image charge repels an adjacent ion. In 
this limit the repulsive interaction between two adjacent 
bound ions is 



Udd[n) 



e^a^it"^ 



/2e. 



(7) 



This repulsion results in the formation of a strongly- 
correlated liquid of ions on the electrode surface, rem- 
iniscent of a two-dimensional Wigner crystal [20|. The 
total dipole energy at a given plate is anSudd(n) , where 
a is a numerical coefficient which depends on the struc- 
ture of the lattice of dipole positions. Thus, the total 
dipole interaction energy can be written as 



Udd = aS[niUdd{ni) + n2Udd{n2)] 



(8) 
4.4. Fur- 



For the minimum energy triangular lattice, a 
ther calculations will use this value. 

We now assemble a full description of the total energy 
U associated with the bound charge densities rii and ri2, 
taking as our reference the case where ni = 712 = 0: 

U = S{ni + n2huim + Udep + UdE + Udd ~ QV. (9) 

Here, —QV represents the work done by the voltage 
source. 

It should be noted that this expression neglects an- 
other potential contribution to the total energy: that of 
the finite chemical potential of ions in the bulk. At zero 
temperature, for example, ions in the bulk are crystal- 
lized into a three-dimensional Wigner crystal with finite 
self-energy uwc ~ —lAe'^N^^^/e per ion [l^. The en- 
ergy U is defined relative to the state where no ions have 
left the bulk, so the attachment of (ni -I- n2)S ions to 
the electrodes should involve an additional energy cost 
of — (ni -I- n2)Suwc- At small iV, this energy provides 
only a small correction to the binding energy Uim and 
can be effectively absorbed into the coefficient 7. 

At equilibrium, the values of Q{V), ni{V), and n2{V) 
are those which minimize U. They can therefore be found 
by the conditions dU/dQ ~ dU/drii = dU/dn2 = 0, the 
latter two of which produce 



5a 



3/2 3 , 



5a 



l'\^ 



47r 


V 


+ 


Q 

eS 


An 


/ 

n-2 
V 


— 


Q 

eS 



: 7+47r 2ni 



zS 



(10) 



7+47r ( 2^2- -^] a^ 
' (11) 



If ions are relatively sparse in the bulk, so that No? <^ 1, 
then the first and final terms of both equations can be 
neglected. Setting Q = suggests that at zero voltage, 
when the net charge of the capacitor is zero, there is still 
a finite concentration no of ions bound to each plate given 
by 



noo 



77Va3 
47r 



<1. 



(12) 



Since the "filling factor" na^ on each metal surface is 
small, bound ions are sufficiently distant that our approx- 
imation of a dipole interaction between them is justified. 
This condition also verifies our assumption that the first 
and final terms of Eqs. (|10l) and PT|) are much smaller 
than unity. 

As the voltage is increased from zero, ions are driven 
away from the anode and attracted to the cathode, so 
that ni decreases and ^2 increases. Eqns. ([TU]) and pT|) 
imply that rii and 71.2 depend linearly on Q: 



no - ni ~ Q/eS, 
n2 - no ~ Q/eS. 



(13) 
(14) 



Subtraction of Eqs. (TTSl) and (fT4|) suggests that the total 
number of adsorbed ions ni -|- n2 — 2no per unit area 
remains almost constant with increasing voltage. This 
result is a direct consequence of the large difference be- 
tween the depletion layer energy Udep and the dipole- 
dipole energy aSnUdd{n)- Indeed, by comparing Eqs. 
([5]) and ([5]), we see that the condition No? <^ 1 implies 
that the energy cost associated with increasing the total 
number of bound ions, and thereby causing the depletion 
layers to swell, is much larger than the dipole-dipole in- 
teraction energy. As a consequence, the total number of 
bound ions remains nearly constant with voltage. The 
electronic charge Q ~ eS{n2 — ni)/2 that passes through 
the voltage source can therefore be thought of as the cor- 
responding movement of image charges from one plate to 
another. 

When enough charge has moved through the voltage 
source that Q — eSriQ, the anode becomes completely 
depleted of bound ions, so that ni = and n2 ~ 2no. 
This corresponds to a particular voltage Vc- At voltages 
V > Vc, the number of bound ions on the non-depleted 
electrode may still increase, but only through the costly 
widening of the depletion layer. 

A relation between charge Q and voltage V can be 
derived by substituting Eqs. (jT^ and (ITU) into Eq. (0), so 
that the total energy U{Q) is written as a function of the 
charge only. A subsequent application of the condition 
dU/dQ — gives 



Vc^- 
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Q 

eS 
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(15) 

By taking the derivative of this expression with respect 
to V, we can derive the capacitance C = dQ/dV. This 



gives 



C = 



15a 



eS 



i^riQa? — Qa^ jeS 



eS 



a^UQa^ + Qa^/eS 



(16) 



In this expression the terms inside the parentheses rep- 
resent the capacitance of the anode and cathode, respec- 
tively, which add Uke capacitors in series. At zero vohage, 
(5 = 0, so that the capacitance becomes 



C(0) 



4 / 47r 
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eS 
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C 



H- 



(17) 

For Na^ < 1 one gets C(0) > Cff/2 because at smah 
voltages charging of the capacitor is Umited only by the 
dipole-dipole repulsion energy. Since uqu^ <C 1, the 
dipole-dipole interaction is weak, so that the resulting 
capacitance can be large. 

At higher voltages, the capacitance of the EDL near 
the positive plate increases strongly as this plate plate 
becomes depleted of ions and the corresponding dipole 
repulsion energy goes to zero. Thus the contribution of 
the positive plate to the total capacitance vanishes when 
Q = enoS. The corresponding voltage Vc can be found 
by substituting Q — en^S into Eq. (|15p : 



(18) 



5a/7iVa3y/^ 



Immediately prior to V = Vc, the capacitance achieves 
its maximum value 



Cri 



1/4 



eS 



= \/2C(0). 



(19) 



15a y^Na?' J a 
The effective thickness <i*jj„ corresponding to Cmax is 



'^™-S(^)'''«"«-4^(^«^)'''- (^«) 



Thus we arrive at a remarkable prediction: the effective 
capacitor thickness can be much smaller than the ion 
diameter a. As an example, an ion-conducting medium 
with mobile sodium atoms of diameter a = 2 A and den- 
sity Na^ = 0.01 can be used to make a capacitor with 
capacitance nearly six times larger than Cr 1'^- As men- 
tioned before, this surprisingly high capacitance is a re- 
sult of the weak dipole-dipole interaction between bound 
ions that comprise the double-layer. Indeed, near the 
capacitance maximum the filling factor on the negative 
plate 7^20^ — 2nQa^ ^ 1, so that it is incorrect to think 
of the EDL as a series of uniformly charged layers. This 
difference represents an important change of paradigm, 
from a mean-field capacitor to a capacitor composed of 
discrete, correlated dipoles. 



At V > Vc, ions can no longer simply be transferred 
from anode to cathode, and the capacitance collapses to 
a much smaller value. The value of this "depleted capac- 
itance" can be found through optimization of the total 
energy U under the condition ni = 0, which yields 



CdepiV) ~ — 



Na^ 



a V MV/Vm + 7) ■ 



(21) 



This expression neglects the weak dipole-dipole interac- 
tion at the non-depleted negative plate. At V/Vim ^ 7 
the capacitance is dominated by the depletion layer next 
to the positive electrode, and therefore it approaches the 
standard value for capacitance of a depletion layer. 

Fig. [3] shows the capacitance and the density of bound 
ions as a function of voltage, as calculated by a numerical 
minimization of the total energy in Eq. (jH]). We have used 
Na^ = 0.1, following the estimate of Ref. [6]. All results 
are within 10% of the approximate analytic expressions 

in Eqs. dUD - (HID, (UH) - dni. 




FIG. 3: The area densities of bound ions (dashed lines, left 
vertical axis) and the capacitance (solid line, right vertical 
axis) as a function of applied voltage for Na'^ = 0.01, as given 
by a numerical minimization of the total energy in Eq. @ 



III. LOW-TEMPERATURE THEORY OF AN 
ASYMMETRIC CAPACITOR 

We may also consider the case of an asymmetric capac- 
itor, where the two electrodes have unequal areas. This 
may be the case, for example, in a coaxial cylindrical 
capacitor where the inner electrode has a smaller radius 
than the outer electrode. For the sake of argument, we as- 
sume that at positive voltages the smaller electrode (area 
^i) is the anode and the larger electrode (area ^2) is the 
cathode. The corresponding Helmholtz capacitance for 
the double layer at each electrode is Ch,i = sSi/2Tra and 
Ch,2 = e5'2/27ra for the small and large electrodes, re- 
spectively. These two capacitances add in series, so that 



by the GCS theory the maximum possible capacitance is 
-1 /I 1 \ -1 



Taking the derivative with respect to V, and using C 
dQ/dV, we get 



1 



1 



Ch,i Ch,: 



— - 
2iia \Si 



1 



(22) 



When 5*1 = 5*2 = 5*, this resuh reduces to ChI'^-, as in 
the previous section. In the hmit S2 ^ S'l, Eq. ((^ 
approaches Ch,i, so that the total capacitance is limited 
by the smaller area. 

Let us now see how our theory of the previous section 
can be generalized to the asymmetric case. The condi- 
tions of neutrality at anode and cathode become 



eNzi — eni — Q/ Si 
eNz2 - en2 = -Q/S2, 



(23) 
(24) 



so that the electrostatic energy of the depletion layers 
can be written as 



^ dep — 



2^ 
3eN 



^^ '- + #;: 



^^'"^"1^ 



(25) 

Similarly, the interaction between dipoles and the deple- 
tion layer becomes 



UdE 



-Anui' 



niSi ni + 



Q_ 

eSi 



+ n2S2 n2 



Q_ 
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(26) 



and the total dipole-dipole interaction energy is 

Udd = a [SiniUdd{ni) + S2n2Uddin2)] ■ (27) 

The total energy is therefore 

U = (51^1 + S2n2hu,^ + Uaep + U^e + Udd - QV. (28) 

As before, the values of Q{V), ni{V), and n2{V), are 
those which minimize U. 

At zero voltage, the conditions dU/drii = dU/dn2 = 
suggest ni = n2 — no, where no is the same as in the 
symmetric case [Eq. P^ ]. At finite V they imply 



no -ni 
712 -"-0 



Q/eSi 
Q/eS'2, 



(29) 
(30) 



so that, as in the symmetric case, the charge Q can be 
thought of as the movement of image charges from one 
electrode to another. Here we have again assumed that 
Na^ < 1. 

Using Eqs. ((29)) and ((30|) . we can rewrite the total en- 
ergy as a function of Q only. By optimizing U{Q) with 
respect to Q we obtain an expression V{Q) from which 
we can define the capacitance. This procedure yields 
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(32) 



In this expression the terms inside the parentheses rep- 
resent the capacitance of the anode and cathode, respec- 
tively, which add like capacitors in series. At zero voltage, 
Q = 0, so that the capacitance becomes 
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(33) 
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When S'l — S2, this expression reduces to the symmetric 
result of Eq. ([T7)) . When Si <^ S2, the capacitance is 
dominated by the smaller area, as expected. In either 
case, C(0) is again much larger than the corresponding 
Helmholtz values Cjj/2 and Ch,i- 

Below we consider the behavior of the capacitance as 
a function of voltage, examining separately the cases of 
positive and negative voltage applied to the small elec- 
trode. 



A. Positive voltage 

When the voltage is increased from zero, the capaci- 
tance increases as ions unbind from the anode and bind 
to the cathode. At a certain critical voltage Vc.i the an- 
ode becomes depleted of bound ions. This occurs when 
7ii = 0, or Q = eSiTiQ, so that by Eq. (PT|) 
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(34) 



At this point, the double-layer capacitance of the anode 
diverges as ions bound to the anode become sparse and 
the corresponding dipole-dipole interaction energy goes 
to zero. The total capacitance is therefore limited only by 
the capacitance of the cathode. Substituting Q — eSino 
into Eq. ^ gives 
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Si 
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cm. 



(35) 
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When the cathode is much larger than the anode, 
S2/S1 3> 1, the maximum capacitance C{Vc^i) — 
C(0)S2/Si. This result implies a tremendous growth 
in the capacitance at positive voltages < ^ < Vc.i- 
For voltages just below the critical value, so that 14,1 — 
y ^ V"c,i, we can examine this growth analytically. At 
S2/S1 :^ 1 we can ignore the term Q/eS'2 in Eq. (1^ 
(31) and we arrive at {noo^ — Qa? /eSi)^^"^ = (K.i — V)/V„n- 



Substituting this result for the first term of the sum in 
Eq. ([5^ (the inverse capacitance of the smaller electrode) 
and ignoring the second term (the inverse capacitance of 
the larger electrode), we find that the capacitance di- 
verges approximately as 



-'"'^ll 



2/3 



V,,r 



Vrl - V 



1/3 



^, (0 < F < K,i) 



(36) 

before being truncated by the finite value of S2 as given 
in Eq. (I55|) . At F = 0, Eq. (pS)) approximately matches 
the capacitance C(0) from Eq. (flT)) . 

At larger voltages V > Vc,i, the anode becomes de- 
pleted of ions and the capacitance collapses to a much 
smaller value. This value can be found through mini- 
mization of the total energy U with respect to Q under 
the condition ni = 0, which gives 



Cdep,l{V) 



Na^ 



eSi 



^TT{V/V,m+"f) a 



{V>Vc,i). (37) 



In other words, the capacitance at large positive voltages 
is dominated by that of the growing depletion layer at 
the anode, as in Eq. (piT) for the symmetric case. 

Fig. m shows the capacitance as a function of voltage 
for different values of S2/S1, as calculated by numerical 
minimization of the total energy in Eq. ([25]) . For V > 0, 
this figure illustrates the analytical results of Eqs. ([M)) 
- (|37p . For y < it shows the predictions of the fol- 
lowing subsection. Notice that at large values of S2/S1 
the divergence in the capacitance near V = Vc,i becomes 
increasingly pronounced. 
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FIG. 4: The capacitance as a function of voltage for different 
ratios S2/S1, as calculated by numerical minimization of the 
total energy, for Na^ = 0.01. Threshold voltages 14, 1 and 
Vc,2 are shown for the dashed line only. 



tion of negative voltage V, where the electrodes change 
roles and the cathode area 5*1 is smaller than the an- 
ode area 82- This change in electrode roles suggests an 
asymmetry in the total capacitance with respect to the 
sign of the voltage. Below we examine the capacitance 
at negative voltages for two limitin g cas es: the moder- 
ately asymmetric case 82/ Si <^ 1/vNa^ and the highly 
asymmetric case S2/S1 » 1/vNa^. 



1. Moderately asymmetric capacitor 

If the two e lectr ode areas are not too different, so that 
S2/S1 <^ 1/vNa^, then the capacitance at negative volt- 
ages can be described using the same procedure as for 
positive voltage, beginning with Eqs. ([5T|) - (P^ . The 
voltage Vc,2 at which the larger electrode becomes de- 
pleted of bound ions can be found by setting Q = —eS2no 
in Ec^. pip , which gives 



K. 



3/2 



Vc,l- 



(38) 



The large electrode therefore becomes depleted at larger 
voltages (in absolute value) than the small electrode (see 
the dashed line in Fig. S]). The condition S2/S1 <C 
1/vNa^ guarantees that |T4,2| < Vim- 

At larger negative voltages V < Vc.2, there is an appar- 
ent collapse in capacitance due to the depletion of bound 
ions at the larger electrode. The capacitance immedi- 
ately prior to this voltage 



J2 



Si 



i^^<°' ^ (I 



3/2 



C(Fc,i)- (39) 



This value is smaller than the maximum capacitance 
C{Vc,i) since the capacitance at 14,2 is limited by the 
capacitance of the smaller electrode while at 14, 1 it is 
limited by the larger electrode. For voltages beyond the 
collapse point, V < 14,2, the total capacitance can be 
found through minimization of U under the condition 
ri2 — 0, which yields 



Cdep,2{V) 



Na^ 



eS2 



47r(7-V^/F„„) a ' 



{V <V,^2)- (40) 



This capacitance is similar to Eqs. (PT|) and ([57|. and is 
determined by the growing depletion layer at the larger 
electrode. 



2. Highly asymmetnc capacitor 



B. Negative voltage 

So far we have been talking about the case of a small 
anode and large cathode. We now consider the situa- 



We now consider the opposite limit of electrode areas, 
when 5*2 and 5*1 are so different that S2/S1 » l/VNa^. 
In this case there is no apparent collapse in the capac- 
itance at any negative voltage — T4m < 1^ < (see the 
dotted line in Fig. g]). Rather, Eqs. ^^ and ([321) imply 
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a continuous increase in capacitance with voltage as in 
Eq. dSni): 



C{V) 



Af2V^f V,m V"" eSi 



3 \5a/ VV^c^i — V J a 

{~y^ra <y < Ka)- (41) 



Notice that at F = —Vim the capacitance becomes ap- 
proximately equal to the Helmholtz value C(— T^m) ~ 

At larger negative voltages V < —Vim our description 
of the total energy loses its validity, since Eq. (pij) im- 
plies a near-complete filling of this electrode by bound 
ions, nio^ ~ 1. Thus the assumption of a dipole-dipole 
interaction between bound ions is no longer accurate. 
At y <C —Vim, the capacitance of the smaller electrode 
should be dominated by the accumulation of complete 
ionic layers at the small ( negat ive ) electrode, as described 
by the theories of Refs. [Il| - [l6j . This is opposite to the re- 
sult of the symmetric case, where the capacitance at large 
negative voltages is determined by the growth of the de- 
pletion layer next to the positive electrode. Instead, the 
large difference between S2 and 5*1 guarantees that the 
capacitance of the smaller electrode is much lower, and 
it therefore determines the total. The corresponding ca- 
pacitance at such large negative voltages is 



C{V) 



— Vim £J1 

AttV a ■ 



{V « -V,m)- 



(42) 



At V = -Vim, this capacitance also approaches the 
Helmholtz value C{—Vim) — Ch,i, so that Eqs. (HT|) and 
(|42| match a.t V ~ -Vm- 

Eqs. (HI]) and ^j indicate that in the limit 5*2/ 5*1 > 
1/y/Na^ the capacitance does not depend on the area of 
the large electrode 5*2. Rather, it is totally determined by 
the physics of the double-layer near the small electrode. 
Therefore, we can view this result as the capacitance of 
a single small electrode (S2 — >■ 00). While the result 
of Eq. (|42p is known [ll| - [l6| . we could not find that of 
Eq. (|41|) in the literature . In the range of its validity, 
the capacitance grows with increasing voltage V from 
C ~ Ch,i at T^ = —Vim and actually diverges at F = 
Vc,i. Thus the effective thickness d* of a single interface 
vanishes at T^ = ^,1! Of course, this result is valid only 
at zero temperature. As shown in next section, finite 
temperature truncates this divergence. 



IV. TEMPERATURE DEPENDENCE OF 
CAPACITANCE 



Wigner crystal. Since the capacitance is limited only 
by the resulting dipole-dipole interaction between ions, 
our prediction for the capacitance is highly dependent 
on the validity of this description. If thermal motion of 
bound ions produces a significant correction to the re- 
pulsive interaction between them, than the capacitance 
will be different than our "zero-temperature" prediction 
of the previous sections. 

In general, when the two electrodes have area S": , 5*2 , 
the total free energy of the system can be written as 



F 






S2n2)lUim 

+ S2f{n2) 



^ dep 

QV, 



UdE 



(43) 



where f{n) represents the free energy per unit area of 
the two-dimensional collection of bound ions that forms 
at a given interface, which includes dipole-dipole repul- 
sion and thermal motion along the metal surface. At 
low temperatures ksT <^ Udd{no), bound ions are crys- 
tallized and f{n) ~ anudd{n), so that Eq. (|43p be- 
comes equal to Eq. (PH)) and we regain the results of 
our low-temperature theory. If the average free energy 
per ion along each metal surface is small enough that 
df/dn <Si Uim, and the voltage is small enough that we 
still have ni,n2 ^ 1/a^, then the equilibrium require- 
ments dF/dni = dF/dn2 = imply 



no- ni = Q/eSi 
112 - no = Q/eS2, 



(44) 
(45) 



as before. 

As a consequence, for sufficiently small voltages \V\ < 
Vci where neither metal surface is depleted of bound ions, 
the sum ni + n2 = 2no remains constant. The condi- 
tion dF/dQ = 0, along with dni/dQ — —1/eSi and 
dn2/dQ = -l/eS2, gives 



eV^f{n2)-f'{m). 



(46) 



Here, the / denotes a derivative with respect to the argu- 
ment. The capacitance dQ/dV is therefore 



C^e" 



f'ini) , r(»2) 
^1 



^2 



(47) 



Generally speaking, Eq. (|47p can be combined with 
Eqs. (|44l) - ([46]) to give an analytic relation for the ca- 
pacitance as a function of voltage at a given temperature: 
C{V, T). This section will focus primarily on the temper- 
ature dependence of the zero- voltage capacitance C(0, T) 
and on the capacitance maximum Cmax{T)- 



In the previous sections, we derive capacitance by min- 
imization of the total electrostatic energy and do not in- 
clude any entropic eflFects associated with the finite ther- 
mal energy of ions. Here we consider the dependence of 
capacitance on temperature. 

Until now, bound ions have been assumed to oc- 
cupy their ground state configuration: a two-dimensional 



A. Zero-voltage capacitance 

According to Eq. (|47p. the capacitance at zero voltage 
can be written 



C(0,T) = 



1 1 
/"(no) V^^^ 



(48) 



so the capacitance is fully determined by the 
(temperature-dependent) free energy of ions along the 
metal surface, f{n). The behavior of this contribution 
can be separated into three regimes of temperature. 



1. Low Temperature 

At very low temperature ksT <^ Uddino), one can 
imagine that each ion in the Wigner crystal undergoes 
small thermal oscillations in the confining potential cre- 
ated by its neighbors. If this potential is expanded to 
second-order in the displacement r from the potential 
energy minimum, then the average squared thermal dis- 
placement is 



ekuT 



IS.Se^a^n^/^ ' 



(49) 



again assuming a triangular lattice of dipoles. The posi- 
tional entropy of a bound ion, relative to an unbound 
state in the bulk, can be estimated as ln((r^)/a'^) — 
ln(l/A''a^). The free energy per unit area f{n) of bound 
ions is therefore 

fin) ~ n [auddin) - fcsTln {N{r^)a)] , (50) 

The corresponding zero- voltage capacitance 

2 ksT 



C(0,r) = C(0,0) 



1 



3a Udd{no) 



(51) 



where C(0, 0) is the zero- voltage capacitance described 
by Eq. HMD. 



2. Intermediate Temperature 

At sufficiently large temperatures that /csT ^ Udd{no) 
but ksT <ti Uim, the crystal-like order of dipoles is de- 
stroyed, and bound ions are better described as a two- 
dimensional ideal gas than as a Wigner crystal. In this 
limit, the free energy per unit area of bound ions f{n) 
can be approximated as that of a two-dimensional ideal 
gas. 
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(52) 



Here we again define the entropy of an ion relative to the 
bulk. 

The description of Eq. ([SS)) assumes that ions are non- 
interacting, so that their free energy is determined purely 
by entropic motion. The effect of relatively weak interac- 
tion between ions can be included by a virial expansion 
of the free energy. 



J^fu + n'kBTBiT). 



(53) 



Here, B{T) is the second virial coefficient, calculated 
from the dipole-dipole interaction energy u[r) between 
bound ions as 



B{T) 
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(55) 
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, the resulting capacitance is 
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(56) 
If the temperature is low enough that fcsT <C 
Uim{Na?'Y/'^, then all other corrections to the capaci- 
tance beyond that of the virial coefficient are parametri- 
cally smaller in No?' . 

One can estimate the transition temperature Ti be- 
tween Eqs. ([STjl and ([56| by equating them, which gives 



Ti « 7.%a^/''udd{no)/ki 



(57) 



As expected, the transition occurs when the thermal en- 
ergy is of the same order as the dipole-dipole interaction 
energy. 



3. High Temperature 

At much larger temperatures kgT ^ u^m, ions no 
longer bind to the metal surface. Since at these tem- 
peratures the change in potential at a given electrode 
is small compared to the thermal energy, the attraction 
of ions to the metal produces only a small perturbation 
in the overall ion density. Ions therefore form a diffuse 
screening layer around each metal surface, with a width 
equal to the Debye-Hiickel screening radius 






(58) 



The resulting capacitance per unit area of the double- 
layer is e/4TTrs, so that the total capacitance of the two 
double-layers in series is 
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47rr, 
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(59) 



The crossover temperature T2 between the Debye-Hiickel 
capacitance Cdh and the intermediate-temperature 
value can be found by equating Eqs. ([55)1 and ([55]) . which 
gives 



T2 = 2"fUim/kl 



(60) 



Fig. [5] shows a schematic depiction of C(0) in all three 
regimes of temperature, plotted for two different values 
of the density Na^. 
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FIG. 5: Schematic depiction of the temperature dependence 
of the zero-voltage capacitance of a parallel-plate capacitor. 
Numbers in parentheses indicate the formula describing each 
portion of the temperature dependence. The dash-dotted line 
corresponds to a higher value of the density Na^ , which pro- 
duces a lower zero-temperature capacitance and a higher ca- 
pacitance at large temperatures. The temperature Ti is indi- 
cated for the solid line only. 



B. Maximum capacitance 

At zero temperature, the capacitance maximum is as- 
sociated with a sharp discontinuity in the C-V curve, as 
shown in Figs. [3] and [D At finite temperature, however, 
entropic effects will tend to smooth out these discontinu- 
ities, thereby shifting both the magnitude and position 
of the capacitance maximum. 

At sufficiently low voltages that neither electrode is de- 
pleted of bound ions and ni,n2 <C 1/a^, Eq. (|47|l implies 
that the capacitance is extremized when 



/"'(«i) = ( I ) r{n,). 



(61) 



This relation can be used to find the capacitance maxi- 
mum Cmax{T) at temperatures T < T2. For high tem- 
peratures T > T2, the capacitance is equal to Cdh and 
is essentially constant in voltage. 

Below we derive the maximum capacitance for the 
symmetric case, S1/S2 — 1, and the highly asymmetric 
case, S1/S2 — >■ 0. 



1. Symmetric capacitor 

When the electrode areas are equal, Si = 82, Eq. (pT 
becomes 



/'"(m) = r'K). 



(62) 



This suggests that at V = 0, where ni = 712 = no, there 
is always either a maximum or a local minimum in the 



capacitance. At low temperatures T <C Ti, V = is a 
local minimum. The maximum can be found by solving 
Eq. (|62p and then substituting the results for ni , 712 into 
Eq. (|i7)) . which gives 



Cmax{T) ~ C„a,(0) 1 + 1.2 



VriQa? 
(T < Ti) (63) 



Here, Cmaa;(0) is the maximum capacitance at zero tem- 
perature, given in Eq. ([TO)) . 

Qualitatively, this result can be explained by consider- 
ing that the capacitance maximum at zero temperature 
is driven by a vanishing dipole-dipole interaction at the 
depleted electrode, which allows the capacitance of that 
electrode to diverge. At finite temperature, the free en- 
ergy of bound ions cannot fall below the thermal energy 
fcsT, so the capacitance of the nearly-depleted electrode 
remains finite. Setting Uddini) = ksT and solving for ni, 
while setting 71,2 = 2noj allows one to derive the result of 
Eq. (|63p to within a numerical coefficient multiplying the 
temperature. Note that as T approaches Ti, the capaci- 
tance maximum approaches C(0,T) as in Eq. (|5ip . 

For larger temperatures T 3> Ti , the capacitance max- 
imum disappears and the function changes concavity 
around V = Q. Thus, the maximum capacitance be- 
comes equal to the zero- voltage capacitance given in Eq. 
([55)1 . A characteristic set of capacitance- voltage curves 
corresponding to this range of temperature is shown in 
Fig.m 



2. Highly asymmetric capacitor 

For the case 6*2 ^ 5*1, our zero-temperature theory in 
the previous section predicts a sharp divergence in the 
capacitance as V approaches Vc.i [Eq. (|36p ]. driven by a 
vanishing dipole-dipole repulsion at the smaller electrode. 
At finite temperature, this divergence is truncated by 
entropic effects, which inhibit the complete depletion of 
bound ions from the electrode surface. According to Eq. 
(j6ip . for S1/S2 — > the maximum is characterized by 



/"'(ni)=0. 



(64) 



By Eq. (|47p . the corresponding capacitance C = 
e^Si/rim). 

At low temperatures T < Ti, the solution of Eq. 
gives a capacitance 

0.32 / U,m \ ^'^ 

a XksT J 



^inax \-^ ) 



^, (r«Tr). 
a 



(65) 



As in the symmetric case, this maximum occurs when 
ni declines sufficiently that Uddini) — ksT. In the limit 
that the temperature T approaches Ti, Cmax — C(0,T), 
and there is no increase in capacitance at positive voltage. 
At intermediate temperatures Ti < T < T2, the ca- 
pacitance C{V) at small voltages \V\ <^ ksT/e is dom- 
inated by the two-dimensional entropy of ions bound to 
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the metal surface 5*1. As a consequence, the capacitance 
increases with negative voltage, where the density of ions 
ni increases and therefore their entropy declines. The 
capacitance continues to rise with negative voltage until 
V ~ —ksT/e, at which point the density of ions is large 
enough that the dipole-dipole repulsion Udd{ni) is compa- 
rable to the ideal gas free energy per ion dfid{ni)/ dni . At 
this point the capacitance achieves its maximum, which 
is again well-described by Eq. (|65|) . Fig. [TOl shows a char- 
acteristic set of C-V curves corresponding to this range 
of temperature. 

At large temperatures T > T2, the maximum ca- 
pacitance becomes similar to the Helmholtz value Ch,i- 
The maximum occurs at large negative voltages V < 
—kBT/e < —Vim, where the applied voltage is strong 
enough to collapse the ionic screening layer to the elec- 
trode surface and form a complete layer. 



Comparison with ion-conducting glass 
experiment 




C /2 

h' 



V (Volts) 

FIG. 6: The capacitance of an ion-conducting glass between 
metal plates. The thin lines with dots show data from Ref. 
[g| from three different samples (reproduced with the authors' 
permission). The heavy soUd line is our theoretical prediction, 
using Na"^ = 0.1, a = 2 A, and 7 = 4. The dotted line shows 
the corresponding Helmholtz value Cii/2. 



Let us return to the case of a symmetric double-sided 
capacitor and compare our theory to the experiments of 
Ref. 6] . Capacitance- voltage characteristics for three dif- 
ferent phosphosilicate glasses are shown in Fig.[6]together 
with our theoretical prediction for the relevant concen- 
tration Na^ = 0.1 and temperature T = 600 K (heavy 
solid line). The theoretical curve is derived by a numer- 
ical minimization of the total free energy, as in Eq. (1431) , 
using e = 10 for the bulk of the glass and e = 2.5 for 
ions bound to the metal surface. This approximation is 
equivalent to using 7 = 4. If one assumes a uniform di- 
electric constant e, or 7 = 1, then the theory predicts 
an even larger capacitance C(0) « 3Ch (see Fig. EJ, but 
it also predicts the capacitance to collapse at a smaller 
voltage than what is seen in experiment. To obtain bet- 
ter agreement with experiment one may need to consider 
the disorder potential acting on mobile ions in the glass, 
but this is outside the scope of the present work. 



V. MONTE CARLO SIMULATION 

In order to test the analytic predictions of the previous 
sections, we present here the results of a simple Monte 
Carlo (MC) simulation that models the behavior of an 
ionic conductor between metal electrodes. As a compu- 
tational convenience, we divide the system depicted in 
Fig. [1] into three slabs and disregard the thick neutral 
middle one, so that more computer time may be devoted 
to the anode and cathode sections in which interesting 
physics is occurring. This separation is in line with the 
above theory, which assumes that the two electrodes are 
separated by a distance much larger than any screening 
length scale, such that there is no interaction between 
the two electrodes' double-layers. 

Each remaining slab is treated as a square prism cell 



X Ly X Lz, where L^ 



Ly and 



with volume fl = L. 
Lz is chosen so that the system is at least twice as thick 
as the depletion layer. The metallic electrode coincides 
with one of the cell's square faces. At F = 0, there is no 
difference between the two partial systems; each contains 
flN mobile positive ions, which are modeled as spheres 
with diameter a — 2A that carry a charge e located at 
their center (the "primitive model"). These mobile ions 
are constrained to move on a cubic lattice with lattice 
constant a, placed so that a/2 is the distance of closest 
approach of an ion to the cell's walls. Mobile ions are 
subject to an excluded volume constraint, so that two of 
them cannot occupy the same site simultaneously. On 
this same lattice, at each of the Mj, = U/a^ total lattice 
sites, there is a small fixed charge —eNa^ which models 
the negative background. The MC program allows for 
Q/e ions to be taken from the anodic cell to the cathodic 
cell without changing the charge of the background. As 
we will see, this movement of Q/e ions is equivalent to 
applying a certain positive voltage V . 

Every charge within a cell forms an electrostatic image 
in the metallic electrode surface (z = 0), i.e. a charge q 
at {x, y, z) has an image charge —q located at {x, y, —z). 
The total electrostatic energy £ of the cell is calculated as 
1/2 times the energy of a system twice as large composed 
of the real charges and their images, so that 



p2 Jl^* 

4e ^, d 






(66) 



Here, qi denotes the charge of particle i, dij denotes the 
distance between particles i and j, and Mt = 2{Mi + Mb) 
is the total number of particles in the system. For real 
mobile ions, qi = e; for the mobile ions' images, qi — — e; 
for the fixed background charges, qi — —eNa^; and for 
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FIG. 7: The average dimensionless ion density pa as a func- 
tion of tlie distance z from the metal surface, as calculated 
from a MC simulation of a 40 x 40 x 40 A^ cell at Q = 0. 
The dotted line shows the concentration which neutralizes 
the negative background. 



the images of the background charges, qi = eNa? . The 
dielectric constant is set to £ = 5 everywhere. 

At the beginning of a MC simulation both the tem- 
perature T and the zero-voltage mobile ion density No?' 
are set. The positive ions are then initialized to random 
non-overlapping coordinates on the lattice and the initial 
energy is calculated from Ec^. (IM]) . After selecting an ion 
at random, the MC program attempts to reposition it to 
a random lattice site within a cubic volume of (4 A)'^ cen- 
tered on the ion's current position. For one in every 100 
attempted moves the MC program expands this volume 
to (20 A)^ in order to overcome the effects of any large, 
local energy barriers. The cell is given periodic bound- 
aries, so that an ion exiting one face of the cell re-enters 
at the opposite face. The total electrostatic energy of 
the system, £", is calculated after each attempted move. 
Moves are then accepted or rejected based on the stan- 
dard Metropolis algorithm. To ensure thermalization, 
8, 000 moves per mobile ion are attempted before any 
simulation data is collected. After thermalization, sim- 
ulations attempt between 10^ and 10^ moves per mobile 
ion in the system, of which ^ 10% — 40% are accepted. 

In Fig. [3 the average ion density as function of the z 
coordinate is shown for a simulation of a 40 x 40 x 40 
A3 cell, with T = 350iC and Na^ = 0.01. The accu- 
mulation of bound ions at the metal surface (z = lA) 
is clear, as is the depletion layer adjacent to it. As pre- 
dicted, the system regains electroneutrality beyond the 
depletion layer, validating our assumption that the two 
double layers of the metallic electrodes can be simulated 
separately. The peaks in ion density at either edge of 
the neutral region likely correspond to over- charging by 
a strongly-correlated liquid of ions in the bulk |23.l23|. 

The primary goal of the simulation is to obtain the full 
system's capacitance, C = dQ/dV. A change in voltage 
in the real system requires charge transfer, Q, from one 
electrode to the other through the voltage source. In 
our simulation we apply a constant voltage 1^ > to 
the pair of cells by adding Q/e ions to one cell (cathodic 



side) and removing Q/e ions from the other cell (anodic 
side) . In this way the electronic charge of each electrode 
is varied; the addition of Q/e ions to a cell corresponds 
to an addition of a charge — Q to the electrode, which 
comes in the form of image charges for the added ions. 
The voltage that corresponds to this movement of charge 
is found through V = dT/dQ, where J^ = Fa + Fc is the 
free energy of the full system, equal to the the sum of 
the anodic and cathodic free energies. In a given cell 
containing Mi positive ions, we used the Widom particle 
insertion method [2J, [23 to obtain the change in either 
the anodic or cathodic cell's free energy AFa^c caused by 
the addition of another positive ion. 



kr,T 



= In {p{x, y, z)aV(exp[-A£(a;, y, z)/{kBT)]}) . 

(67) 

Here, A£{x, y, z) is the change in electrostatic energy due 
to a probe charge placed at {x,y,z) and p{x,y,z) is the 
mean density of real positive ions at (x, y, z). The angle 
brackets denote a time average and the ratio of the quan- 
tities inside the natural log is independent of the position 
{x,y,z). At the beginning of the simulation, 20 lattice 
sites are selected as Widom insertion sites at which the 
quantities p{x, y, z) and exp[— A£(x, y, z)/{kBT)] are cal- 
culated after every attempted move. The corresponding 
values of AF^ j, obtained from each insertion site are av- 
eraged in order to give a final value. 

We are free to set the free energy of the system aX Q — 
to zero, so that J-"(0) = 0. Simulating the anode and 
cathode for Q — 1,2,3..., while employing the Widom 
insertion method, allowed us to determine the system's 
free energy as a function of Q, 



T{Q) = J-(Q - 1) + AF,(0) - AFaiQ). 



(68) 



This equation must be used iteratively to find T{Q — 1) 
starting from F{Q) = 0. For the symmetric case {Si = 
S2), this process is easily extended to negative voltages 
by taking Q < 0; positive ions are then attracted to the 
anode and repelled from the cathode. In the asymmetric 
case {SI < S2), AFc{Q) = AFa{-QSi/S2), so we need 
only to simulate the cathodic cell at positive and negative 
Q in order to calculate AFa{Q) and F{Q). 

Taking the discrete derivative of these data points gives 
J-{Q) — J-{Q — e) = V{Q — e/2). Another derivative gives 
the capacitance of the system as a function of V, 



C{V{Q)) = 



dQ e 

dV " {V{Q + e/2) - V{Q - e/2)) ' 



(69) 



where V{Q) = [V{Q + e/2) + V{Q - e/2)]/2. 

Fig. [8] shows the results of the simulation for a system 
with 5*1 = 5*2 = 40 X 40 k^, L^ = 20A, and Na^ = 
0.01. Capacitance as a function of voltage is shown for 
three values of the temperature, along with the analytic 
predictions explained in section ITVl These temperatures 
fall within the range Ti ^ T < T2, and as predicted, 
the maximum capacitance occurs at zero voltage while 
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are also shown. Results are plotted as a function of ab- 
solute temperature as well as dimensionless temperature 

T* = kBT/{e^/ea). 



4- 



I 
O 3 



2- 



o 

<M 1 



I 
O 3 



o 

(M 



- b) □ 



0.05 



T* 



0.1 



0.15 



' 






"^ \ 


Eq.(56) 

— ■ — numerical solution 
D MC simulation 


- 


Na^ = 0.01 


= 


— 



1,000 



Na-" = 0.1 



2,000 



3,000 




1,000 2,000 

I (K) 



3,000 



FIG. 8: Capacitance as a function of voltage for a symmetric 
parallel-plate capacitor at various temperatures Ti <T < T2, 
using Na'^ = 0.01, e = 5, and 7 = 1. a) Results from a Monte 
Carlo simulation of the OCP model. Error bars for the MC 
data are smaller than the symbols, b) Analytic predictions, 
as explained in section ITVl 



the capacitance collapse is smeared over a voltage range 
proportional to fc^T/e. 

In order to quantify the finite-size effects of our simu- 
lation cell, we examined the capacitance at zero voltage, 
C(0,T), obtained from three "slab-shaped" simulation 
volumes of size L x L x L/2, with L = 40, 60, and 80 A. 
For No?" = 0.01, C(0,r) was seen to scale linearly with 
1/L at all values of the temperature that we examined 
(T = 350, 600, and 1200 K). In each case, the value of 
C(0,r) obtained by extrapolation to infinite system size 
was within 16% of the value of C(0, T) corresponding to 
i = 40 A. This difference was within the uncertainty of 
our simulation for all temperatures except 1200 K. We 
also checked that there was no dependence of the capaci- 
tance on the aspect ratio of our simulation cell by exam- 
ining C(0, T) in three cubic cells with side length L — 20, 
40, and 60 A. The resulting value of C(0, T) again scaled 
linearly with 1/L and, within uncertainty, the extrapo- 
lated values of C(0, T) agreed with those found in the 
slab geometry. These results allow us to conclude that 
the 40 X 40 X 20 A simulation cell provides a good approx- 
imation of an infinite system. All MC results presented 
below correspond to this choice. 

The temperatures explored by our MC simulation fall 
in the intermediate temperature range Ti ^ T < T2, and 
so the capacitance should be described by Eq. ([5S|) . In- 
deed, as shown in Fig. [9l the temperature dependence of 
C(0,T) obtained from simulations is in good agreement 
with the analytical prediction of Eq. (|56|) . Results from a 
numerical minimization of the total free energy [Eq. (1431) ] 



FIG. 9: Capacitance at zero voltage, C(0, T), as a function of 
temperature for a symmetric system Si = ^2 at two different 
values of the density Na"*. Squares represent the MC data 
using S = 40 X 40A'^ , while the prediction of Eq. (|56|) is shown 
by the thick black curve. Results are plotted as a function of 
absolute temperature (lower axis) as well as dimensionless 
temperature T* — ksT /{e^ /ea) (upper axis). The error bars 



for the MC data are smaller than the symbols, a) Na 
b) Na^ = 0.1. 



0.01. 



Unfortunately, the low-temperature predictions of our 
theory, corresponding to T ^ Ti, could be not be ex- 
amined directly since these temperatures correspond to 
an extremely low acceptance rate in our MC simulation. 
Nonetheless, we can get an idea of the zero-temperature 
capacitance by examining the behavior of the total elec- 
trostatic energy £. Since at zero temperature the to- 
tal free energy becomes equal to £, the capacitance C 
approaches {(P£/dQ'^)~^ at low temperatures. Examin- 
ing {(P £ / dQ^)^^ as a function of temperature and ex- 
trapolating to T = allows us to make a rough esti- 
mate of the zero-temperature capacitance C(0,0). For 
Na^ = 0.01, the resuh is 2C(0,0)/Cif w 7, which is 
significantly higher than the theoretical prediction of 
2C(0,0)/Ci/ = 4.5 given by Eq. ^. This enhanced 
capacitance may be the result of screening of the dipole 
interaction by mobile ions in the bulk, which suppresses 
the interaction of distant dipoles and therefore reduces 
the effective value of a. At larger ion density, the dis- 
crepancy between our low-temperature theory and the 
projected zero-temperature capacitance from simulation 
becomes even more pronounced. For Na? = 0.03, 0.1, 0.3 
we estimate 2C(0, 0)7(7// = 6.5 ± 1, as compared to 
2C(0,0)/Ci/ = 3.5,2.5,2.0 given by Eq. ^. These 
results are consistent with the interpretation based on 
screening of dipole-dipole repulsion by ions of the bulk. 
Indeed, at larger ion densities the bulk becomes more ef- 
fective at screening because it is separated from the metal 
surface by a thinner depletion layer. 
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As a rudimentary test of the extent to which the dipole 
interaction is screened by bulk ions, we performed a MC 
simulation of a 40 x 40 x 20 A simulation cell in which 
one of the ions was fixed to the center of the metal sur- 
face (a; = y = 0, z = a/ 2). The time-averaged density of 
ions p{x, y, z) was then recorded at every lattice site in 
the simulation cell, from which the mean electric poten- 
tial (j){x, y,z = a/2) at the metal surface could be recon- 
structed. In order to isolate the contribution of bulk ions 
to the potential from that of strongly-correlated ions on 
the metal surface, the simulation was run with naS ions 
removed from the system. In this way we simulated the 
anodic side of a capacitor at V = T4, where the metal 
surface is depleted of bound ions. At a density corre- 
sponding to No?' — 0.01 and a temperature T = 350 K, 
we found that the potential cj) surrounding the ion was 
well described by (f){r, z — a/2) = ea^ /2er^ •exp(— r/r^fc), 
where r = \/ x^ + j/^ is the azimuthal distance from the 
bound ion and Vsh is a length scale which characterizes 
the range of the screened dipole potential. We found 
Tg ss 4.5a, which is smaller than the average distance 

— 1/2 

riQ « 6a between bound ions at zero voltage. If this 
expression e0(n~^/^) is substituted for the dipole-dipole 
energy M^d, then the resulting prediction for zero- voltage 
capacitance at No?" = 0.01 is significantly enhanced: 
2C(0)/Cff == 12 as compared to 4.5 from Eq. ^^. In 
reality, the observed capacitance from MC simulations is 
between these two values, 2C{Q)/Ch ~ 7, and this dis- 
crepancy may be the result of a non-additive response of 
bulk charges to dipoles at the surface (non-linear screen- 
ing). 

Finally, we also considered the simulation of a highly 
asymmetric capacitor, where 5*2 is infinite and 51 = 40 x 

40 A . We again examine the case A^a^^ = 0.01, T ~ 
600K, and use Lz = 20 A. The results are shown in Fig. 
[TOl along with the analytic predictions of section IIVI for 
different values of the temperature. 



VI. LOW- VOLTAGE CAPACITANCE PEAK IN 
ASYMMETRIC IONIC LIQUIDS 

Up to this point we have discussed the OCP model 
primarily as it applies to ion-conducting glasses at mod- 
erately high temperatures where only the smallest ion 
(cation) is mobile. Qualitative agreement of our predic- 
tions with experimental data shows that the OCP model 
is a reasonable zero-order approximation. 

In this section we discuss other possible applications 
of the OCP model, as mentioned in the introduction. 
One application which immediately comes to mind is 
to super-ionic crystals, where only the smallest positive 
ionic species (such as Na"*" or Li+) is mobile. In this case 
there is no reason for a frozen disorder, so that in super- 
ionic crystals the OCP model should work even better 
than in ion-conducting glasses. 

Perhaps more interesting is the application to ionic 
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FIG. 10: The capacitance of a highly asymmetric capacitor, 
S2/S1 — !■ cxD, where only surface 1 determines the capacitance, 
a) Results from a MC simulation. Error bars are smaller than 
the symbols, b) Analytic predictions for various tempera- 
tures, as explained in section IIVI Note that the capacitance 
diverges at T — >■ 0. 



liquids, which have recently attracted considerable at- 
tention [1^ I26l - l28l |. In ionic liquids, both positive and 
negative ions are mobile. In order to get spontaneous po- 
larization near the electrodes at zero voltage, and thereby 
obtain the low-voltage peak in capacitance predicted by 
the OCP model, one should consider a strongly asym- 
metric ionic liquid. 

We begin by considering ionic liquids composed of 
monovalent cations, which we model as rigid spheres 
with diameter a, and much larger monovalent anions, 
modeled by rigid spheres with diameter A ^ a. For 
example, one may have in mind the ionic liquid made 
of Na"*" cations and large non-coordinating anions such 
as the "BARF" anion ([B[3,5 - (Ci^3)2C6i/3]4]~) ^■ 
In such a liquid anions fill most of the space and form 
a weakly-compressible negative background. The maxi- 
mum energy of attraction of the anion to its image in the 
metal surface, e'^/2eA, is much smaller than the corre- 
sponding image attraction e^/2ea for a cation. Because 
of their small size, cations easily move between anions. 
The cations are strongly attracted to the metal plates and 
therefore rearrange themselves to form the EDL. Thus, 
we may assume that anions form an analog of the nega- 
tive background in the OCP model. The maximum den- 
sity of such a background is approximately 1/A'^. Thus, 
we can expect that the capacitance C(V, T) is similar 
to what we predict in the OPC model, if for N we use 
N ^ 1/^'^, or in other words Na^ — > [a/A)^. For exam- 
ple, if A/ a = 100"'^''^ = 4.6 we should get a capacitance 
C(0,r) similar to the case Na^ = 0.01 studied above for 



15 



the OCP model. 

In order to verify these predictions, we ran MC sim- 
ulations of the primitive hard sphere model of an ionic 
liquid between equal-sized metallic electrodes. In these 
simulations, monovalent cations and anions are given a 
diameter a and A, respectively, and placed with a partic- 
ular volume density iV in a square prism simulation cell 
with volume fl = L x L x L/2. The electrode surface is 
again chosen to coincide with the z = plane, so that an 
ion of charge q = ±e whose center is at position {x, y, z) 
has an image charge —q at {x, y, —z). The voltage of the 
electrode is varied by changing the number of cations 
Mc in the system, as in simulations of the OCP model, 
while the number of anions Ma — VlN remains constant. 
The corresponding electronic charge in the electrode is 
Q = e{Ma — Mc) and the capacitance dQ/dV can be 
determined from the resulting voltage. A very similar 
simulation method was used by previous authors 13 Ct |3l| 
to examine the capacitance of ionic liquids, but the effect 
of asymmetric ion size was not explored. 

Since in this case anions are mobile, unlike in the OCP 
model, a change in voltage should correspond to a chang- 
ing number of anions as well as cations in the vicin- 
ity of the metal surface. Therefore, one may object to 
our method of modifying the charge of the electrode by 
changing only the number of cations in the system. For 
example, one may imagine inducing a charge Q = —2e 
in the electrode by democratically adding one cation and 
removing one anion from the simulation cell rather than 
by adding two cations. The difference between these two 
methods, however, is only an infinitesimal change in the 
bulk ion densities of the simulation cell; the physics of the 
metal interface is not affected. To ensure that the capac- 
itance in our simulation is independent of the method of 
charge transfer, we repeated all of our simulations using 
both a method where the charge Q is modified by chang- 
ing only the number of anions Ma and a method where 
Ma and Mb are changed simultaneously by equal and op- 
posite amounts. No noticeable change was observed to 
any of the results presented below. 

The microscopic rules of the simulation are identical 
to those of the OCP model, except that ions are not 
constrained to move on a lattice and there is no fixed 
negative background. The energy of a particular config- 
uration of ions is also identically calculated, with the ex- 
ception that the hard-core repulsion between ions should 
be added explicitly to Eq. (|55)) . That is. 






(70) 



where Mt = 2 (Ma -I- Mc) is the total number of charges 
in the system (ions plus images), and the two-particle 
interaction energy uidij) is 



.^. . loo, d,,<(D.,+D,)l2 

"■ ''" Iq.g.M,, d,,, >(A + ^,)/2' ^ ' 



Here, Di denotes the diameter of ion i; Di — a for cations 
and Di = A for anions. 

In addition to the method of Eq. ([M)) for calculat- 
ing capacitance, where the voltage is inferred from the 
change in free energy of the system, for these simula- 
tions we used a method where the voltage is measured 
directly for a given value of Q, so that determination 
of the free energy is unnecessary for calculating capaci- 
tance. The voltage is measured by defining a "measure- 
ment volume" near the back of the simulation cell — 
occupying the range — L/4 < x < L/4, — L/4 < y < L/4, 
i/4 < z < 3L/8 — inside of which the electric poten- 
tial is measured. After performing thermalization of the 
initial random configuration (50,000 moves per mobile 
ion), the total electric potential (/)(x, y, z) is measured at 
500 equally-spaced points within the measurement vol- 
ume after every 3(Ma -|- Mc) attempted moves. These 
measured values of potential are then averaged both tem- 
porally and spatially to produce a value for the elec- 
tric potential 4>{Q) of the electrode relative to the bulk. 
The corresponding voltage between the two electrodes is 
V ~ 4>{Q) — 4>{—Q), and the capacitance is determined 
from the discrete derivative AQ/AT/. The results pro- 
duced by this second method were compared with those 
produced by the method of the previous section for four 
different sets of simulation parameters, and the results 
were indiscriminable. Below we present results from only 
the second, more time effective method. 

Fig. [TT] shows the resulting capacitance for ion liquids 
with asymmetry A/a = 4 and A/a = 2 at various tem- 
peratures. Here the temperature is presented in dimen- 
sionless units T* — ksT / { e^ / ea) in order to facilitate 
comparisons with literature [30|-|32| . Our simulations use 
the same values for the cation size and dielectric constant 
as in the OCP case, a = 2 A and e = 5, so that the tem- 
perature scale e^ /eaks ~ 16700 K and the range of data 
0.02 < T* < 0.18 corresponds to 350 K < T < 3000 K. 
For a more typical value of the cation diameter a ~ 8 A, 
this range corresponds to 100 K < T < 750 K. The size 
of the simulation cell for the A/ a = 4 case was L = 80 A, 
and in the case A/ a = 2 we used L = 40 A. The dimen- 
sionless ion density p^ = Ma{a^ -I- J^)/^ was p\ ~ 0.5. 
An examination of finite size effects, as in the previous 
section, suggests that our results for capacitance are ac- 
curate to within 18%. As predicted above, our numerical 
results for ionic liquids with A/ a — 4 and 2 are close to 
the results for the OCP model with Na^ = 0.01 and 0.1, 
respectively. 

While Fig. [TT] shows the capacitance of the two- 
electrode system, the capacitance of a single interface 
can also be easily determined from our MC simulations 
by looking at the derivative d4>/dQ. As an example, Fig. 
1121 shows this capacitance as a function of dimensionless 
voltage V/Vim for the case A/a = 2 at a temperature 
T* = 0.036 and density pi = 0.5. The asymmetry in 
capacitance with voltage is similar to what we observed 
in the OCP model (Fig.[Tni). 

The results of Fig. [T^ can be compared to those of a 



16 




b) 



O 3 



o 

(M 



A/a = 2 



0.05 



0.1 

T* 



0.15 



0.2 



FIG. 11: The capacitance of a primitive model ionic liquid 
between metallic electrodes as a function of temperature. The 
density of ions in each case is pi = 0.5. Error bars are smaller 
than the symbols, a) A/a = 4. b) A/a — 2. Compare these 
results to those of the OCP model in Fig. [9] 
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FIG. 12: The capacitance of a single metal/ionic liquid inter- 
face as a function of voltage for an ionic liquid with A/a = 2. 
The temperature is T* — 0.036 and the density is pi = 0.5. 
Error bars are smaller than the symbols. Compare these re- 
sults to those of the OCP model in Fig. [TO] 



recent study [33| which performed a molecular dynamics 
simulation of an ionic liquid with size ratio A/ a — 2 
near a uniformly charged interface. In this study, the 
capacitance was found to be similarly asymmetric with 
voltage, with the larger capacitance resulting when the 
voltage has the same sign as the larger ion. However, 
since the study considered only the response of the ionic 
liquid to a uniformly-charged plane (essentially treating 
the metal electrode as a charged insulator), it arrived at 
capacitances C < Ch at all voltages. It is also worth 
emphasizing that over the range of voltage in Fig. [T^] 
we are still far from complete coverage of the electrode 



surface. Indeed, at V/Vim = 0.5 the area coverage of 
the surface is only 35% and 0.2% by anions and cations, 
respectively, while at V/Vim = —0.5 it is 3% and 28%. 
Thus, in our case the collapse in capacitance is not caused 
by the building of a second layer of ions, as proposed by 
Ref. [11. 

Thus far we have focused on the effects of asymmet- 
ric ion size, but we note that there is another way to 
make a strongly asymmetric ion liquid. Namely, cations 
and anions may have the same radius a but different 
absolute values of charge. One can imagine, for exam- 
ple, that cations are multivalent and have charge +Ze 
while anions have charge — e. In this case cations are 
much more strongly attracted to their —Ze images, so 
that together they again create a dipole layer on the sur- 
face of the metal. Because there are Z anions per one 
cation, anions form a thicker negatively-charged layer 
centered farther from the metal surface than the cations. 
This anion layer is analogous to the depletion layer of in 
the OCP theory, with a dimensionless concentration of 
cations Na^ ~ l/(-^ + !)■ In order to estimate the ca- 
pacitance, we can use the results of the OCP model with 
Na^ ^ 1/{Z+1). 

This prediction can be checked by our MC methods by 
simulating an ionic liquid with trivalent cations (Z = 3) 
and a neutralizing concentration of monovalent anions, 
both with the same diameter a. We consider the case 
where the dimensionless temperature, normalized to the 
larger charge of the cation, is T* = ksT /{Z'^e^ /ea) — 
0.11 and the density is p^ = 0.5. The resulting capac- 
itance is shown in Fig. [13] as a function of the dimen- 
sionless voltage V/Vim, where Vim = Ze/2ea. The C-V 
curve is again very similar to that of the OCP model, 
with a maximum C /Ch > 1 and a smaller capacitance 
at positive voltage, where Z-ions are depleted from the 
electrode surface. 
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FIG. 13: The capacitance of a single metal/ionic liquid inter- 
face as a function of voltage for an ionic liquid with trivalent 
cations and monovalent anions with equal diameter. The di- 
mensionless temperature T* — ksT /{Z'^e^ /ea) = 0.11 and 
the density is pi = 0.5. Error bars are smaller than the sym- 
bols. Compare these results to those of the OCP model in 
Fig.[TOl 



Even closer imitation of the OCP model can likely be 
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obtained by combining a large charge asymmetry Z ^ 1 
with a large asymmetry of ion size A/ a 3> 1. 



VII. AQUEOUS SOLUTION OF A Z:l SALT 

Another system for which the OCP model gives a zero- 
order approximation is that of an aqueous solution of a 
salt with multivalent cations with charge Ze ( "Z-ions" ) 
and monovalent anions, for example LaCla. In such a so- 
lution Z-ions are strongly attracted to their —Ze images 
in the metal electrode and so a number of them form 
compact dipoles at the interface. On the other hand, at- 
traction of anions to their images is much weaker, so that 
at room temperature they do not form dipoles but rather 
stay in the solution and effectively form a negative OCP 
model background with charge density — eZiV, where N 
is the concentration of salt. If we again use the theory of 
Sec. In] to balance the depletion layer energy with the en- 
ergy of Z-ions condensed on the metal surface, we arrive 
at Eq. (fT2|) for the surface concentration no of Z-ions. 
Remarkably, both no and the capacitance at T = T^ = 
do not depend on Z (one can guess this from the fact that 
the corresponding results in section |TI] do not depend on 
the elementary charge e). 

For LaCla we can examine the case where A^ = 0.5 M 
when the salt is totally dissociated. This concentration 
corresponds to Na^ ~ 0.1, if for the diameter of the hy- 
drated La+^ ion we use a — 6 A. Thus, it is tempting to 
apply to this case the above finite temperature calcula- 
tions of capacitance for Na^ = 0.1. 

One may worry, however, that the depletion layer, 
where the concentration of La"'"'^ ions vanishes, is not 
uniformly charged by Cl~ ions with their average den- 
sity — 3eiV. This may happen because the positive po- 
tential (ps near the surface of the metal, which devel- 
ops to balance the image attraction of Z-ions to the 
metal, is much larger than fcsT/e and therefore results 
in the exponential growth of Cl^ concentration near the 
metal surface. For the surface potential we get Ze4>s = 
Z'^e'^ /2ea — kBTh\{na? /Na^) , where the second term 
comes from the entropy that a Z-ion loses at the surface 
in comparison with the bulk solution. For LaCla, using 
a = 6 A, e = 80, r = 300 K, and No? = 0.1, we get that 
Z'^e^/2ea = b.SksT and \n{noa'^/Na^) < 1, and there- 
fore, e(j}s is equal to only LlksT . This allows us to ignore 
(as a zero-order approximation) the non- uniformity of the 
concentration of Cl^ in the depletion layer. Then we can 
use the temperature-dependent results we obtained for 
the OCP model in Section IIVI Applying these results 
requires only the scaling of the temperature in units of 
Z^e^/ea. For LaCls, using e = 80 and a = 6 A, we 
find that the temperature unit Z^e^/ea is 0.37 of that 
for glass (where Z = 1, e = 5, and a — 2 K). Thus, for 
7V~0.5MatT = 300 K we arrive at the same ratio 
C{0)/Ch — 2.5 for the capacitance of a single interface 
as for the OCP model with Na^ = 0.1 at T = 800 K 

(Fig.Eb). 



This prediction can be checked by running a MC simu- 
lation identical to the one described in the previous sec- 
tion. Fig. [TJ] shows the resulting capacitance per unit 
area C/S of a single interface as a function of voltage for 
the salt concentrations N — 0.5 M and A^ = 1.5 M, using 
the temperature T = 300 K and the estimated hydrated 
diameters a = 6 A and A = 4 A for La"*""^ and Cl~, re- 
spectively. The simulation cell is given a size L = 100 A. 
For N — 0.5 M and 1.5 M, the maximum capacitance 
is larger than the Helmholtz value by 2.0 and 2.3 times, 
respectively. The data presented in Fig. [U corresponds 
to the range of electrode charge | AQj < ZeNil./2, where 
ZeNil represents the total cation charge in the simu- 
lation volume. Restricting our simulation to this range 
ensures that the bulk ion concentration is not changed 
significantly by the addition/removal of Z-ions to the cell 
that is associated with finite voltage. 
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FIG. 14: The capacitance per unit area C/S of an aqueous 
solution of LaCla at a metal electrode, as determined by our 
MC simulation using the primitive hard sphere model. The 
La^"^ cation has a hydrated diameter a — 6 A, while the Cl~ 
anion has a diameter ^ = 4 A. The temperature is T = 300 
K and the uniform dielectric constant e — 80. Error bars 
are smaller than the symbol size. The dotted line shows the 
Helmholtz value Ch/S ~ 236 /iF/cm^. Compare these results 
to those of the OCP model in Fig. \W\ 



The results of Fig. [14] are qualitatively similar to those 
of the OCP model given in Fig. |TU] The capacitance 
is asymmetric with respect to voltage, acquiring a larger 
value at negative voltages where Z-ions accumulate at the 
metal surface and bind strongly to their image charges. 
At positive voltages, Z-ions are depleted from the elec- 
trode surface and the capacitance collapses. At large 
positive voltages V ^ kBT/e « 26 mV, negative an- 
ions become strongly bound to the electrode by the ap- 
plied voltage and they approach complete filling of an 
ionic layer. In this limit the capacitance approaches 
e/2TrA = Ch ■ a/ A. 
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VIII. CONCLUSIONS 

In this paper we have presented a theory to explain 
how the capacitance of the metal/ionic conductor inter- 
face can be significantly larger than the Helmholtz ca- 
pacitance Ch- In other words, we have shown how the 
apparent thickness of the double layer d* can be smaller 
than the ion radius. This surprising conclusion is ob- 
tained by abandoning the mean-field approximation and 
considering instead the behavior of discrete charges next 
to the metal surface. While mean field theories cannot 
explain how d* can be smaller than the physical separa- 
tion between the electrode and its countercharge, we have 
shown that very large capacitance is a natural result for 
an EDL composed of discrete, correlated ions. We have 
worked within the approximation of a "one-component 
plasma" model, where only one species of ion is mobile, 
and described its behavior over the full range of temper- 
ature. We have further argued that our results can be 
easily extended to strongly asymmetric ionic liquids. A 
simple Monte Carlo simulation confirms our analytical 
predictions at realistic temperatures. At very low tem- 
peratures, the EDL capacitance is limited only by the 
weak repulsion between ion-image dipoles at the metal 
surface, which for a single interface produces a sharp ca- 
pacitance peak that diverges as T"^/'^. 

Qualitatively, our theory explains all the main features 
of the experiment in Ref. [y] . The authors of Ref. Q re- 
late their observations to theories of so-called "pseudo- 



capacitance" , a term used by Conway and coworkers for 
the rare cases of anomalously large EDL capacitance (see 
[2] and references therein). Pseudo-capacitance is said to 
result from specific adsorption of cations to the metal 
surface, where the cations are neutralized. In this sense 
Conway's theory is similar to ours. However, his theory of 
pseudo-capacitance does not explain what happens with 
the negative charge of excess anions, which remains in the 
bulk and which in our theory plays a pivotal role. The 
existing theory of pseudo-capacitance also does not ex- 
plicitly specify the form of the repulsion between bound 
ions, and therefore does not arrive at a closed result. We 
take care to address both of these points in the present 
(OCP) model, and we arrive at definite predictions for 
capacitance. Thus, one may consider our theory to be an 
improved theory of pseudo-capacitance, if by this term 
one understands a capacitance larger than the Helmholtz 
value. We emphasize, however, that our theory does not 
assume any Faradaic effects, so that our result is in fact 
a standard capacitance and the prefix "pseudo-" is un- 
necessary. 
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